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Abstract 

The distribution functions of photons in squeezed and correlated light for one-mode and 
multimode cases are obtained based on the method of integrals of motion. Correlation 
coefficient and squeezing parameter are calculated. The possibility to generate squeezed light 
using nonstationary Casimir effect is discussed. Quantum parametric Josephson junction is 
proposed as quantum vacuum generator of electrical vibrations. 


1 Introduction 


The aim of this work is to discuss integrals of the motion and uncertainty relations and to obtain 
the distribution function of photons in squeezed and correlated light for one-mode and multimode 
cases. The distribution function of photons in squeezed light for one-mode fields was discussed 
by Schleich and Wheeler [1], by Agarwal and Adam [2], and by Chaturvedi and Srinivasan [3]. 
The photon distribution function for squeezed and correlated light [4] and [5] was discussed by 
Dodonov, Klimov and Man’ko [6]. This distribution function depends not only on the squeezing 
parameter, but also on the correlation parameter connected with Schrodinger uncertainty relation 
[7] as well, 


SqSp > 



( 1 ) 


where the parameter r is the correlation coefficient of the position and momentum 


r = (SqSp) 1 1 — - (g)(p)| ■ (2) 

The states with nonzero parameter r we call the correlated states. In the section below we’ll 
consider the problem how to find the states which minimize the Schrodinger uncertainty relation. 
For such states instead of the Schrodinger inequality we have the equality 


SqSp — 



(3) 


These states describe squeezed and correlated light. We will demonstrate in the next section 
how these states are naturally created for quantum parametric oscillator. The case of the photon 
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distribution function for the two- mode squeezed light was considered by Caves, Zhu, Milburn and 
Schleich [8]. Multidimensional generalization of the expression for the distribution of photons in 
squeezed light in terms of Hermite polynomials of several variables may be reformulated. We 
derive this expression on the bases of the result obtained in [17], [9] and [10] for a nonstationary 
parametric multidimensional oscillator. 

The squeezing phenomenon in quantum optics is closely related to the oscillator models de- 
scribed by relativistic wave equations for elementary particles with mass spectrum. These rela- 
tivistic models have been studied by Yukawa [11], by Markov [12], by Ginzburg and Man’ko [13], 
and by Kim and Noz [14]. As shown in [14], the Lorentz boost applied to relativistic oscillator 
gives the squeezing whose mathematics is identical to that of the squeezing in quantum optics. 
The statistical properties of such squeezed relativistic oscillators have been studied by Kim and 
Wigner [15]. 

To obtain the photon distribution function we will consider the nonslationary multidimensional 
oscillator. We shall discuss first the one-mode case in Sec. 2. 


2 One- mode Light 

The Hamiltonian for one-mode light is given by the formula 

H = huj{a)a + ^). (4) 

This mode of the electromagnetic field in a resonator may be described by the model of the 
mechanical oscillator with the Hamiltonian 


p 2 muj 2 q 2 

H = 2 m + 2 

In this case the annihilation and creation operators with boson commutation relations 

P \ 


where 
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(7) 

( 8 ) 


connect both Hamiltonians and forms, together with the identity operator, the basis of the 
Heisenberg- Weyl algebra. In coordinate representation the complete set of coherent states | a) 
satisfying the equation 

a | a) = a | a), (9) 

where a is any complex number, is given by the formula 


{q | a) = 7r *1 2 exp 
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The dispersions of the positions 8q and the momentum 8p do not depend on the parameter a, and 
are given by the relation 

vT 

P 0 
v/2' 

For the coherent states the product of these dispersions minimizes the Heisenberg inequality 


8q = 

8p = 


(11) 

( 12 ) 


kip = 2 * 


(13) 


The time evolution of the coherent state | a, t) may be obtained by simple replacement of the 
parameter a in the formula (10) by the term aexp( — iut) and the phase of the wave function. We 
have 

— iujt 

(q I a,t) = {q | a exp (— zu><)) exp ( - ). (14) 

The correlation coefficient of the position and momentum is equal to zero for arbitrary coherent 
state. It is also equal to zero for stationary Fock state | n,t) satisfying the eigenvalue equation 


a} a | n, /} = n | n,i), n = 0, 1, 2, ... 

This state has the following wave function in the coordinate representation 


(q | n, t) — 'K */ 22 ?(n!) *H n (j)ex p 


q - . 1 * 
-W- l “ t{n+ 2 ] 


(15) 


(16) 


The photon distribution function W n {a ) for the coherent state | a, t) is determined by the overlap 
integral 

|(n,t|a,0| 2 = ^ n (a) (17) 

and coincides with the Poisson distribution function 

I p n 

W n {oi) = - - exp (-| a | 2 ). 


nl 


The mode has the following time-dependent integral of the motion 

A(t) — exp ( iujt)a . 


(18) 


(19) 


We now discuss how the influence of the dependence of the oscillator frequency fi(£) on time 
will change the photon distribution function and the dispersions of the conjugate variables q and 
p . The Hamiltonian of the mechanical parametric oscillator depends on time and has the form 


m-f 

2m 2 

This system has the linear integral of the motion [16] 




( 20 ) 


(21) 
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Here u> = fl(0), and the complex function t(t) satisfies the equation of classical oscillator motion 

e + tt 2 (i)e = 0. (22) 

The integrals of motion A(t) and A^{t) satisfy the boson commutation relation 

[A(t)J\t)] = 1, (23) 

if the Wronskian for the equation (22) is given by the equality 

ce* — e* e = 2 iu>. (24) 

The initial condition for the function e(t) may be taken as follows 

e(0) = 1, e(0 ) = ioj. (25) 

If the frequency of the oscillator is constant, the function e — exp (iujt) and the formula (21) gives 
the integral of the motion (19). The normalized state | 0 ,t) satisfying the Schrodinger equation 
and the relation 

A(t)\0,t)=0 (26) 

has the following wave function in the coordinate representation 

(?|0,0 = 7r ‘*( /c )"^ ex p(s?)- ( 27 ) 

The state | a,t) which is the eigenstate of the integral of the motion A{t) given by formula (21) 

A(t) | a,t) = a | a,t) (28) 

has the following wave function in the coordinate representation 

{q \a,t) = (q\ 0, t) exp a T + • (29) 

Here a is an arbitrary complex number and 

(I3,t | a,t) = exp ( 30 ) 

The wave function (29) satisfies the Schrodinger equation. The Fock states of the parametric 
mode | n, t) satisfy the eigenvalue equation 

A*(t)A(t) | n,t) = n | n,t), n = 0,1,2,... (31) 

The solutions to this equation have the following form in the coordinate representation 

(,|„, ( ) = ^(,|0 >i )(|A„( r i? rl ). (32) 
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Since the state | a,t) is the generating state for the Fock states | m,<) 

| a | 2 ~ a m |m,<) 

\a,t) = exp ( —) 2_s 7=7— » 

2 m=0 yjm\ 

the transition probability from the initial state | n) may be calculated 

WT = rn > n. 

772! \ 2 / 


( 33 ) 


(34) 


Here the transition probability W$ is the probability to be in the ground state 

Wq = 2(| e | 2 + u> _2 | e | 2 + 2) _1/2 . (35) 

For n > m the formula (34) must be changed to 

W7 = ^fp^(K)) 2 . ( 36 ) 

The numbers n and m in the formulae (34) and (35) are either both even or both odd. If one of 
these numbers is even and another number is odd the transition probability between such states 
is equal to zero 

W% +1 = W % + 1 = 0, k,p = 0, 1,2, ... (37) 

The formulae (34) and (36) describe the photon distribution function for the one-mode electromag- 
netic field in a resonator either with moving walls or with media with time-dependent refraction 
index. Thus, we conclude that the squeezing parameters of the parametric oscillator 



(^)’Vm, 

(38) 


/ 2 \ 1/2 

(39) 

s p = 

: LJ * p=|e| 


are connected with the photon distribution function by the ratio (35) which may be rewritten in 
the form 

Wq = 2(S ? 2 + S 2 p + 2)~ 1/2 . (40) 

In the case of vacuum light S p = S q = 1 and the vacuum-vacuum transition probability Wfj is 
equal to unity. 

Another photon distribution function corresponds to the excitation of light state which may 
be described by the model of the forced mechanical oscillator with Hamiltonian 

m = (41) 

This oscillator has the integral of the motion [16] 

A(t) = exp (iut)a + 6(t), (42) 


133 


M = 


(43) 


^WiJo /( T ) eX P(*WT)rfT 


If the initial state of a forced oscillator is the coherent state, the squeezing parameters S q and 
S p are time-independent. They are equal to unity. The photon distribution function is described 
by a Poisson distribution. Thus, if the initial state is the vacuum state the Poisson distribution 
has the form 

| 8 | 2m 


i / i r 1 2 \ 


(44) 

m! exp(| 8 j ) 

The physical meaning of the parameter | 8 | 2 (43) which determines the integral of the motion (42) 
is just the mean photon number after the excitation of the vacuum state by the external linear 
force. The photon distribution function W™ in the case when the initial state was the state | n) 
with n photons is described by the function 


W? = 


n 


H 8 


|2(m— n) 


m! ex 


p(im 


/-:■"(! < i 2 )] 2 . 


(45) 


Here the function L™ is the Laguerre polynomial 

Now consider a general situation when the frequency of an oscillator depends on time and an 
external force is present. The Hamiltonian of the mechanical oscillator model looks like 


H(t) = ~P 2 + \u 2 {t)q 2 - f{t)q- 

We have taken m = uj — h ~ 1. The linear integral of motion A ( / ) is equal in this case to 

,4(i) — 1 / ( / ) fi T T £(i), 

«(0 = ^( e (0 - «( 0 )i 

v(t) = - (e(<) + ic(t )) , 


= — l -jK J /(T)c(r)dr. 


v/2 

The normalized eigenstate x /’ 0 (^,t) of the integral of motion (47) has the form 

|2 \/2aq ' a(<5e* + 6*e) c* 2 e * 


xl’ a (q,t) = </’ 0 (<M) ex P 
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where 


V’o(<7, t) = 


exp 
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(46) 

(47) 

(48) 

(49) 

(50) 

(51) 

(52) 


The squeezing parameters S q and S p for the states (51) are described by the formulae for the 
unforced parametric oscillator (38), (39). The correlation coefficient r is given by the expression 


r=l«|-' [(«) ! -l]'T 


(53) 
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The Fock states which are the eigenstates of the integral of motion A^(t)A(t) are of the form 




( 54 ) 


The photon distribution function for the electromagnetic field created due to the nonstationary 
Casimir effect is expressed in terms of Hermite polynomials of two variables 


where 


W™ 2 


«• , 

= X2 = 6-—, 


(55) 


(56) 


and the matrix R has the elements 


i? = r 1 


T] -1 

-1 -if 


The parameters ( and 77 are given by the relation 

c(t) = <e“ - ,e-“. 


(57) 


The photon distribution function (55) has oscillatory behavior due to the oscillatory behavior of 
the Hermite polynomial of two variables. 

The last photon distribution function describes the influence of the nonstationary Casimir 
effect on the initially thermal equilibrium state 


/5(0) = Z 1 exp 


0{a ] a + 


1. 


(58) 


Z " 1 = 2sinh(/?/2). (59) 

The distribution of photons in the light mode is expressed by the density matrix diagonal elements 


where 


and 


R = D 


p — = -Uw(*,,* 2 ), 

Poo n\ 

_ j f 7/*u(l — Z 2 ) —Z 

—z uu*(l — z 2 ) , 


(60) 


8*uz + 8v* 


|2’ 


X 2 = X C. 


_ U I Z* “ I V 

The photon distribution function has the following deformed Planck distribution form 

1 . ,2^ + 1 


(n) = 
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(61) 
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If there is no external force the parameter <5 = 0. The deformed Planck distribution has the form 

(n)--^- + |H 2 coth(/?/2). (62) 

e p — 1 

The squeezing parameters S q and S p depend on temperature 

S, =| f | (coth(/3/2)) l,! , (63) 

S„=|<l(coth(/?/2)) 1/2 . (64) 

Thus the parametric excitation of the oscillator may produce the squeezing phenomenon when 

Sq < 1, (65) 

or 

S p < 1. (66) 

But the higher the temperature the more difficult to obtain the squeezing. 

3 Polymode Squeezed Light 

We will consider the photon distribution function for polymode squeezed and correlated light 
using the model of nonstationary parametric multidimensional quantum oscillator with N degrees 
of freedom. Its Hamiltonian may be written in the form 

ff = iqB(()q + <?(<)q. (67) 

where the vector 

q = = (p,x) (68) 

contains N momentum projection operators and N position projection operators. The 2Nx2N 
- matrix B(t) and 2N - vector C(t) are time-dependent parameters of the system. The model 
corresponds to N light modes in the resonator. The interaction of these modes depends on time 
either due to the motion of the resonator walls or due to the time-dependence of the media 
refraction index. The system must demonstrate the properties of nonstationary Casimir effect for 
N t— ► oo . The oscillator has 2N - vector 

Q (t) = A(t)q + A(t), (69) 

which is the linear integral of motion if the 2Nx2N-symplectic matrix A satisfies the classical 
equation of motion 

A = (70) 


where the 2Nx2N-matrix E has the form 


— # 

The 2N - vector A (2) obeys the equation 

A =■ A (t)E<?(<). 


(72) 


The solution (A, A) of equations (70) and (72) describe the classical trajectories of multidimen 
sional oscillators and may be considered as an element of the inhomogeneous real symplectic group 
ISp(2N,R). The initial conditions for these equations are 

A(0) = 1, A(0) = 0. (73) 

The propagator of the system has the form 

1 f l 

G(x 2 , Xi,t)= 7 exp { --(x 2 AJ 1 A 4 x 2 

■^/det(— 2 z 7 tA 3 ) 2 

— 2x 2 A 3 'xi T xjAjAg ’xi T 2x 2 A 3 S 2 


+2x 1 (^i - \xK%) + S 2 AiAj 1 ^ - 2 fifyr) j , (74) 

with (h = 1). Here the matrices A,, i = 1,2, 3, 4 are NxN blocks of the matrix A 

a =G: t) (75 > 

— # — # — * 

and N-vectors £i,< 5 2 , are the components of the vector A 

A = (£,£)• (76) 

The Hamiltonian of the system may be rewritten in terms of the boson annihilation and creation 
operators a = (aj, ..., a^), a* = (aj, ajy), in the form 

H = i(a,& ')D(t) ( * t ) +/a + fa f . (77) 

If we introduce the 2N - vector 

A(0=(? t ) (78) 

this vector is connected with the 2N - vector q = (p,x) by the relation 

A - Vq, (79) 

where the 2Nx2N - matrix V has the form 

v =M-i !)• (80) 

Then the matrix B(t) in (67) is connected with the matrix D(t) in (77) by the relation 

D{t) = V'BVl (81) 
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The integrals of motion of the multimode nonstationary system are given by the formula 

A (‘) = (iO =Ai(i) (&t) +f(< )- 

The 2N - vector T(t) = (7, 7*) is connected with the vector A 

m = v& (t). 

The 2Nx2N - matrix Ai(t) with NxN - blocks £ and 77 

A ( t \ - ( ^ ’/(O ^ 

l() vrw at)) 

is connected with the matrix A (t) by the relation 

A i(t) = VA{t)V'. 

The propagator of the system in coherent state representation has the Gaussian form 

G(a*J,t) = 


ex P 1 T]d* + a*£ d*( 1 7 +^ 7 ^ */? 


\/det £{t) 


+/3(7* - rj*£ '7) + ^7 I 7~ / l{ T )l( T ) dl 
z Jo 


( 82 ) 


(83) 


(84) 


(85) 


(86) 


The photon distribution function for multimode case may be obtained by expansion of the prop- 
agator (86) into a series with respect to the parameters a*. We have the distribution of photons 
in squeezed and correlated light 


I <m | /M) | = WM = | G( OJ) 


TTl\ ! . . . 771 yv 
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!!-'-! L- 


f' [fi- 


ll 


(87) 


Thus, the distribution function for N - mode system is described by the Hermite polynomial of N 
variables. For squeezed and correlated light the behavior of the function W^(/9) is very oscillatory 
as well as for one- and two- mode cases. The partial cases for two- mode light may be obtained if 
one uses the formulae for Hermite polynomials of two variables found in [5]. 


4 Nonstationary Casimir Effect 

Now let us discuss some possible applications. One of the possible methods to generate squeezed 
light is to use the nonstationary Casimir effect when moving resonator walls produce continuous 
time-dependent reconstruction of the electromagnetic vacuum state. The work against the Casimir 
forces produces two effects. The first effect is the generation of photons from the vacuum. Thus, 
the resonator with mechanically trembling walls is a quantum vacuum generator of electromagnetic 
radiation. The second effect is the squeezing of quantized modes in the resonator due to parametric 
change of vacuum energy. Both effects exist simultaneously. Thus, the plates in the Casimir effect 
may be moved by external mechanical forces. The refraction index of the media may vary with 
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time, the geometrical dimensions of the resonator may be influenced by external mechanical forces. 
In all these cases the vacuum state energy must be changed. This means that the vacuum state 
is continuously reconstructed. For each of the time moments the state is no longer the vacuum 
state due to the change of the parameters. If we have a system with photons or other quanta in 
a box with changing volume, it corresponds to the process with the creation of quanta, e.g. the 
photons. So, due to the work against the Casimir forces one form of energy may be converted 
into the other form. Thus for waving neutral plates (due to mechanical external forces) between 
the plates the photons must be created and this means that the mechanical energy from external 
sources is converted into electromagnetic energy of photons. It is interesting that this effect must 
create the quanta of all other fields existing in nature. Thus, due to Casimir forces we can have 
the generation of photons in a parametric resonator which may be called a quantum vacuum 
generator. 

It is possible to discuss another reduction of nonstationary Casimir effect using the Josephson 
junction. If there is no external voltage in the Josephson junction but its parameters are time- 
dependent, the vibrations of current and voltage will be excited in it. This suggestion [17] is based 
on the analogy of the Josephson junction and a conventional resonant circuit (quantized resonant 
circuit). In classical resonant circuits it is impossible to excite electrical vibrations without external 
sources of voltage. But for a quantum resonator circuit due to Casimir nonstationary forces, it 
is possible to transform mechanical energy which may be the reason for the change of the circuit 
parameters into electrical energy of current vibrations. If this idea is realized it will be a quantum 
vacuum generator of electrical vibrations. The current and voltage in this case play the role of 
conjugate quantum observables and in parametric Josephson junctions they may be squeezed. 
Thus, the quantum noise in Josephson junctions may be reduced for current. In this case the 
voltage will have larger noise. The squeezed and correlated states of Josephson junctions may be 
also excited by changing its parameters with time. 
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